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1  ABSTRACT 

Standard  numerical  techniques  for  removing  aircraft  motion 
and  discontinuities  from  airborne  laser  profiles  are  adapted  to  a 
desk-top  computer.  Because  such  a  computer  is  much  slower  than  a 
larger  machine,  analog  active  filters  replace  the  numerical 
(Hanming)  filters  normally  used.  The  computer  programs  are  in 
BASIC,  and  listings  are  provided  in  an  Appendix.  Several  exam¬ 
ples  of  data  editing  procedures  are  given. 
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BACKGROUND 


Airborne  laser  profiling  remains  by  far  the  most  rapid  method  of  obtaining  some 
measure  of  the  "roughness"  of  a  field  of  ice.  This  paper  will  deliberately  not  deal 
with  the  applications  of  laser  profiles,  but  will  describe  a  specific  method  of 
obtaining  "clean”  profiles  from  raw  data.  The  reader  is  referred  to  Ketchum  (1971) 
for  a  general  description  of  the  profiling  system,  a  Spectra-Physics  Geodolite  3A. 

Ketchum  (1971)  Introduced  airborne  laser  profiles  to  the  Arctic  remote  sensing 
community,  and  Hlbler  (1972)  devised  a  computer-based  scheme  for  removing  unwanted 
"aircraft  motion”  from  the  raw  data.  This  scheme  has  been  used  and  modified  in 
subsequent  applications  (Holyer  et  al.,  1977;  Lohanick,  1979),  but  not  greatly 
changed. 

This  paper  will  describe  a  specific  data  reduction  system  written  at  NORDA, 
centered  around  Hewlett-Packard  hardware  and  firmware.  It  is  based  on  the  Hibler 
three-step  filter  process,  but  is  accomplished  on  a  mini-computer  dedicated  to  this 
task,  and  with  some  steps  which  greatly  reduce  overall  data  reduction  time. 

EVALUATION  OF  RAW  PROFILE  DATA 


Profile  data  are  recorded  in  the  field  on  FM  analog  tapes.  A  time  code  channel 
is  provided  for  comparison  of  the  data  with  annotations  on  various  aircraft  logs, 
and  for  correlation  with  other  airborne  sensors  which  also  record  time.  These 
analog  tapes  must  first  be  converted  to  digital  form  (digitized)  for  handling  in  the 
computer.  The  initially  digitized  data  from  the  analog  field  tapes  will  hereafter 
be  referred  to  as  the  "raw"  data. 

For  those  unfamiliar  with  the  "look"  of  raw  laser  profile  data.  Figure  1  is 
provided  and  annotated  to  point  out  relevant  features.  The  vertical  exaggeration  in 
this  plot  is  about  300:1. 

SAMPLING  RATE 

The  sampling  rate  (number  of  points  digitized  per  second)  for  this  type  of  data 
has  historically  given  an  actual  spacing  of  about  1  m  (see  Hibler,  1972;  Holyer  et 
al.,  1977)  on  the  terrain.  Hibler 's  (1972)  scheme  leaves  intact  all  data  that  have 
a  wavelength  less  than  about  127  m,  although  the  power  spectrum  of  the  topside  of 
sea  ice  has  never  specifically  been  discussed  in  the  short  wavelength  region  (see 
Hibler  and  LeSchack,  1972). 

AIRCRAFT  MOTION 

Long-wavelength  undulations  in  the  data  are  due  to  the  motion  of  the  aircraft 
carrying  the  laser.  As  the  aircraft  slowly  rises  and  falls  in  altitude,  and  makes 
turns,  the  apparent  altitude  measured  by  the  profiler  varies.  This  result  of 
aircraft  instabilities  must  be  removed  if  the  actual  terrain  profile  is  to  be 
recovered. 

DISCONTINUITIES 

Discontinuities  in  the  apparent  profile,  such  as  at  points  A  and  B  in  Figure  1, 
can  be  caused  either  deliberately  or  accidentally  by  the  profiling  system.  Causes 
for  these  jumps  in  the  apparent  profile  are  discussed  in  Holyer  et  al.  (1977). 

Figure  2  shows  the  results  of  not  removing  these  discontinuities.  Nonexistent 
terrain  features  appear  in  final  profile,  because  the  computer  recognizes  only  long- 
wavelength  variations  as  aircraft  motion. 
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Hibler's  1972  scheme  for  aircraft  motion  removal  is  adequate,  and  lends  itself 
to  numerical  or  analog  processes,  as  will  be  discussed  below.  Removal  of  disconti¬ 
nuities  is  essential  to  the  data  reduction.  But,  as  pointed  out  by  Holyer  et  al. 
(1977),  no  satisfactory  method  has  yet  been  devised  for  reliably  eliminating  these 
discontinuities  automatically  (i.e.,  without  the  intervention  of  a  trained  analyst). 
Therefore,  some  sort  of  interactive  system  must  be  devised  to  allow  a  trained 
operator  to  choose  sections  of  data  visually,  as  well  as  to  remove  undesirable 
discontinuities.  The  alternative  to  interactive  graphics  is  a  laborious  process 
involving  hand  annotation  of  computer-generated  plots  and  print  listings  (Lohanlck, 
1979). 


Figure  3  is  included  to  show  a  data  record  for  which  reliable  automatic 
algorithms  for  locating  discontinuities  would  be  difficult  to  construct.  The  data 
were  recorded  with  four  times  the  vertical  sensitivity  of  those  in  Figure  1.  The 
circled  area  is  referred  to  in  a  later  section. 

RELEVANT  CAPABILITIES  OF  THE  COMPUTER  HARDWARE 


All  elements  of  the  hardware  system  used  in  the  present  application  are 
sketched  in  Figure  4.  The  computer  is  a  Hewlett-Packard  (HP)  9845B  with  186  kilo¬ 
bytes  of  internal  memory,  a  CRT  for  printing,  data  display  and  graphics,  an  internal 
thermal  printer  and  two  data  cartridge  (tape)  drives.  The  9845B  has  four  interface 
busses  for  peripherals. 

CRT  graphics  is  done  on  a  560  x  455  dot  matrix.  Any  plot  generated  on  the  CRT 
can  be  printed  on  the  internal  printer  with  one  command.  It  is  also  possible  to 
store  any  plot  on  a  mass-storage  device  such  as  a  flexible  disk  or  tape  cartridge. 

The  HP2240A  Measurement  and  Control  Processor  contains  both  an  A/D  converter 
and  a  D/A  converter,  whose  uses  will  be  described  later.  Both  converters  have  a 
resolution  of  5  mV  in  a  range  of  +10,000  mV  (approximate  range).  The  2240A  is 
directly  compatible  with  the  984 5B.  The  data  transfer  rate  for  A/D  (12-bit) 
conversion  is  about  100  Hz.  The  HP9885  Flexible  Disk  drives  use  seven-inch  disks 
(floppies)  with  a  storage  capacity  of  about  500  K  bytea. 

Two  double  Rockland  active  filters  (Models  432  and  452-01)  are  used  between  the 
2240A  D/A  output  and  A/D  input  to  accomplish  aircraft  motion  removal.  They  replace 
previously  used  numerical  filters,  which  are  far  too  slow  in  this  application.  A 
30000-data  point  record, ' which  would  require  60  hours  of  time  on  this  computer  with 
the  numerical  filters,  requires  about  1.5  hours  with  active  filters. 


THE  DATA  REDUCTION  PROCESS 
DIGITIZING  (PROGRAM  "A/D") 

Because  of  hardware  and  firmware  constraints,  a  maximum  of  90,000  regularly-spaced 
integer  data  points  (16  bits)  can  be  taken  from  an  analog  field  tape  at  one  time. 
Since  the  maximum  sampling  rate  of  the  system  is  100  12-bit  samples/second,  this 
amounts  to  15  minutes  of  data,  or  about  50  nautical  miles  (at  200  knots  aircraft 
speed). 

A  data  tape  is  chosen  and  placed  on  the  Ampex  FR-1300  Tape  Recorder.  A 
suitable  section  of  the  tape  is  located  by  reading  the  time-code  channel  and 
comparing  with  aircraft  logs. 


2 


The  computer  program  used  is  named  ”A/D“  and  its  action  is  diagrammed  in  Figure 
3  (all  programs  are  in  BASIC,  and  complete  listings  are  provided  in  Appendix  A,  with 
internal  documentation).  The  program  enables  the  HP2240A  to  read  the  analog  tape, 
digitize  the  analog  voltages,  and  make  them  available  for  reading  by  the  HP9845B 
computer.  For  maximum  data  transfer  rate  (100  Hz)  the  2240A  constantly  fills  its 
output  buffers  with  data,  while  the  9845B  asynchronously  empties  them  into  internal 
memory.  In  this  way,  time-consuming  synchronization  is  avoided. 

Once  the  90000  (or  fewer)  data  points  are  read,  digitizing  ceases,  and  the  data 
is  immediately  stored  on  a  floppy.  Because  firmware  restricts  any  file  to  32767 
elements  (2^-1),  a  file  length  of  30000  is  chosen.  Thus,  a  full  run  will  fill 
three  files,  each  containing  five  minutes  of  data. 

During  program  execution,  the  operator  presses  one  key  to  begin  or  end  the 
digitizing,  and  names  the  files  which  will  store  the  raw  profile  data. 

REFORMATTING  AND  PLOTTING  (PROGRAM  "PLOT") 

Some  reformatting  of  data  is  done  to  take  advantage  of  special  dense  data  file 
types  which  allow  storage  of  eight  full  files  (40  minutes  of  data)  on  each  floppy. 
The  data  in  one  file  are  then  plotted  in  small  scale  on  the  CRT  and  stored  on  disk. 

A  hard  copy  of  the  CRT  plot  is  obtained  on  the  thermal  printer  (see  Fig.  6),  and  is 
used  for  some  annotation  and  time-scaling.  The  stored  plot  is  used  for  editing  the 
data  as  described  below. 

The  small  scale  plot  is  the  visual  counterpart  of  a  raw  data  file,  and  both  are 
saved  permanently. 

Program  "Plot"  requires  the  operator  to  name  the  file  which  will  store  the 

plot. 

EDITING  (PROGRAM  "EDIT") 

Program  "Edit”,  takes  about  85%  of  total  operator  time,  since  it  requires  him 
to  make  judgments  as  to  the  quality  and  character  of  the  data  which  has  been 
digitized  and  plotted.  He  must  also  choose  one  of  several  editing  procedures  to 
remove  undesirable  discontinuities.  At  this  point,  the  operator  must  be  a  trained 
analyst  who  has  seen  and  been  guided  through,  much  raw  profile  data  and  its 
peculiarities.  The  operator's  expertise  will,  of  course,  have  a  direct  impact  on 
the  quality  of  the  final  product.  Some  examples  given  below  will  help  to  point  this 
up. 


To  begin,  the  operator  names  the  raw  data  file  and  its  plot  (file).  The 
computer  reads  these  into  memory  and  graphics  memory,  respectively,  and  presents  the 
snail-scale  plot  (Fig.  6,  again)  on  the  CRT.  The  operator  uses  a  small  cursor 
(operated  from  the  keyboard)  to  move  to  features  he  wishes  to  change  or  to  observe 
more  closely.  He  presses  one  key,  and  the  CRT  provides  a  large-scale  plot  of  60 
data  points  in  the  immediate  vicinity  of  the  current  point  of  interest  (Fig.  7). 

This  large-scale  plot  will  be  called  an  edit  window.  The  vertical  exaggeration  in 
the  edit  window  is  about  2:1. 

Five  different  editing  procedures  are  available.  Two  of  these  require  the 
operator  to  mark  the  beginning  and  end  of  the  data  record  (BOR,  EOR),  which  are  not 
defaulted.  The  three  editing  procedures  which  affect  the  data  are: 

1.  Straight-line  tie  (see  Fig.  8a).  This  procedure  allows  the  operator  to 
remove  a  sharp  discontinuity  and  replace  it  with  a  horizontal  straight  line.  It  is 
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appropriate  when  the  discontinuity  is  a  "180°  phase  shift"  (an  end-of-scale  jump 
back  to  center  scale,  provided  by  the  laser  hardware  to  avoid  having  to  record 
voltages  outside  a  selected  range,  usually  +1*5  volts),  or  a  spike  which  may  be 
caused  by  spurious  effects  during  recording  or  digitizing. 

The  operator  places  the  cursor  to  the  immediate  left  of  the  discontinuity, 
presses  a  key  to  enter  its  position,  moves  the  cursor  to  the  present  position  of  the 

end  of  the  feature,  presses  the  key,  moves  the  cursor  to  the  desired  position  of 

this  second  point  (only  up  and  down  motion  is  permitted),  and  enters  it.  The  CRT 
shows  the  result  for  approval  or  disapproval  (Fig.  8b). 

2.  Straight-line  slope  (see  Fig.  9a).  The  operator  marks  all  the  points  as 
above.  The  computer  leaves  the  character  of  the  data  in  the  interval  between  points 
1  and  2  intact,  but  changes  the  slope  (Fig.  9b). 

3.  Polnt-by-point  fix  (see  Fig.  10).  Allows  the  operator  to  move  each  point  a 
desired  distance  up  or  down.  This  procedure  is  appropriate  when  there  are  several 
discontinuities  in  the  edit  window,  and  the  data  are  considered  reliable  by  the 
operator.  Several  steps  in  this  process  are  shown  in  Figures  10b  through  lOd. 

The  product  of  this  program  is  a  tape  file  of  all  edit  windows  for  a 

particular  record  of  raw  data  (which  may  be  a  complete  raw  data  file  or  portion  of 

it,  delineated  by  the  BOR  and  EOR  mentioned  above).  This  program  is  the  only  step 
requiring  great  skill  on  the  part  of  the  operator. 

REMOVING  THE  DISCONTINUITIES  (PROGRAM  "DISCON") 

The  operator  provides  the  raw  data  file  disk  and  the  correct  edit  window  file, 
and  the  computer  automatically  removes  all  discontinuities  as  prescribed  in  the 
editing  procedure  above.  The  operator  provides  a  name  for  the  updated  file,  which 
is  stored  again  on  a  floppy.  The  original  file  is  not  destroyed.  Edit  window  files 
are  generally  not  saved,  since  they  are  intermediate  products  and  are  not  needed 
beyond  this  stage. 

AIRCRAFT  MOTION  REMOVAL  (PROGRAM  "8-IN-l") 

The  operator  provides  the  program  with  the  names  of  the  source  file(s)  (the 
one(s)  created  in  Program  Discon  above)  and  the  result  file(s).  The  program  is 
written  to  handle  up  to  eight  pairs  of  files  at  a  time. 

Aircraft  motion  removal  has  previously  been  done  with  numerical  Hamming  filters 
(Hibler,  1972;  Holyer  et  al.,  1977;  Ketchum,  1971;  and  Lohanick,  1979).  The  approx¬ 
imate  number  of  multiplications  to  be  done  by  one  pass  of  a  Hamming  filter  is  the 
number  of  points  in  the  record  times  the  number  of  filter  weights.  For  the  low  pass 
filter  providing  the  aircraft  motion  envelope  (see  below),  the  number  of  filter 
weights  for  a  sampling  rate  of  100  points/second  is  about  250,  leading  to  over  six 
million  calculations  for  a  30000-point  record.  On  the  9845B  this  would  take  many 
hours.  In  fact,  for  all  three  filters,  estimated  total  run  time  for  a  five-minute 
record  is  about  60  hours,  so  we  must  use  external  active  filters. 

Data  flow  through  a  particular  filter  is  diagrammed  in  Figure  11.  Data  are 
first  read  into  an  array  in  computer  memory.  Then  one  point  at  a  time  is  read  into 
the  D/A  converter  and  fed  to  the  appropriate  active  filter.  The  output  of  the 
filter  (connected  to  the  A/D  converter)  is  "simultaneously"  read  back  into  computer 
memory.  Data  are  filtered  in  pseudo-real-time,  since  one  complete  cycle  of  D/A  and 
A/D  conversion  is  about  six  times  as  long  as  the  actual  data  sampling  interval. 
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Thus,  a  30000-point  record  (5  minutes  o£  data)  takes  about  3C  minutes  to  filter. 
Filter  cutoff  frequencies  must  be  set  to  about  one-sixth  of  the  values  calculated 
for  the  numerical  filters.  The  entire  aircraft  motion  removal  scheme  is  explained 
in  Hibler  (1972). 

Since  eight  files  (40  minutes  of  data)  can  be  handled  by  the  program  at  one 
time,  a  complete  run  takes  about  12  hours,  and  is  run  unattended  overnight. 

Operator  time  is  a  total  of  about  15  minutes  in  setting  up.  Result  files,  which  are 
the  final  terrain  profiles,  are  stored  on  disk. 

One  product  of  Program  "8-in-l"  is  shown  in  Figure  12.  This  scaled  plot 
(horizontal  dotted  lines  are  at  2  m  intervals  above  the  zero  line)  has  a  great 
amount  of  vertical  exaggeration  (about  1000:1).  The  operator  uses  this  plot  to  look 
for  any  unusual  features,  such  as  ridges  over  10  m  high  or  large  spikes  below  the 
zero  line.  Causes  for  these  features  oust  be  found.  If  the  cause  is  faulty 
editing,  then  editing  must  be  repeated,  but  only  for  this  particular  error. 

The  profile  is  now  in  its  final  form  on  disc,  and  is  ready  to  be  used  in  any 
application  of  sea-ice  terrain  features  such  as  roughness,  ridge  counts,  ridge 
height  distribution,  or  power  spectral  density. 

A  summary  of  all  steps  taken  in  the  data  reduction  is  shown  in  Figure  13,  for  a 
segment  about  2  km  long. 

REFERENCES 

Hibler,  W.  D.  Ill  (1972).  Removal  of  Aircraft  Altitude  Variations  from  Laser 
Profiles  of  the  Arctic  Ice  Pack,  J.  Geophys.  Res.,  77(30)  p.  7190-7195. 

Hibler,  W.  D.  Ill  and  L.  A.  LeSchack  (1972).  Power  Spectrum  Analysis  of 
Undersea  ar.d  Surface  Sea  Ice  Profiles.  J.  Glaciol.  11,  63,  p.  345-356. 

Holyer,  I.  J.  J.,  P.  Wadhams  and  R.  T.  Lowry  (1977).  An  Interactive  Graphics 
System  for  the  Reduction  of  Airborne  Laser  Profiles  of  Sea  Ice.  Scott  Polar 
Research  Institute  Tech.  Rept.  77-1,  Cambridge,  England. 

Ketchum,  R.  D.,  Jr.  (1971).  Airborne  Laser  Profiles  of  the  Arctic  Pack  Ice, 
Remote  Sensing  of  Environment  2,  p.  41-52. 

Lohanlck,  A.  W.  (1979).  Airborne  Laser  Sea  Ice  Profiles  Near  a  Drifting  Camp, 
April  1977,  NORDA  Technical  Note  49. 


5 


Rough  terroin 


Figure  1 . 


Laser  sea-ice  profile  showing  stages  of  discontinuity  and  aircraft 
motion  removal.  Vertical  scale  exaggeration  about  300:1 . 


Original  Profile 


Final  Profile  with  Discontinuities  not  Removed 


5  meter 
"ridge" 


3.  Laser  sea-ice  profile  with  many  discontinuities.  Circled  area 
is  one  of  several  in  which  two  or  three  discontinuities  occur  in 
rapid  succession. 


igure  4.  Physical  diagram  of  computer  used  in  laser  profile  reduction 
(disc  storage  not  shown) 


Data  Acquisition 


Figure  6.  Small-scale  CRT  plot  of  5-minute  laser  profile  record.  Shown 
actual  size. 
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Figure  7. 


Large-scale  plot  of  60  points  in  the  vicinity  of  a  profile 
discontinuity.  Shown  actual  size.  Vertical  exaggeration 


2.25:1. 
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(a)  Plot  of  discontinuity  as  presented  on  CRT 


orrected  portion 


(b)  Plot  of  removed  discontinuity 


Correct 


" — 4-  _  _ _ 

Straight-line  segment 


Figure  8.  Straight-line  editing  of  a  profile  discontinuity  on  the  CRT. 

Shown  actual  size.  Vertical  exaggeration  of  profile  2.25:1 


(a)  Plot  of  discontinuity  as  presented  on  CRT 


STRAIGHT-LINE  SLOPE 


(b)  Plot  of  removed  discontinuity 


Figure  9.  Slope  removal  of  a  profile  discontinuity  on  the  CRT  shown 
actual  size.  Vertical  exaggeration  of  profile  2.25:1 . 


(c)  Second  discontinuity  removed 


Figure  10.  Point-by-point  removal  of  discontinuities  from  sea-ice 
profile.  Each  frame  actually  fills  CRT.  Vertical  scale 
exaggeration  10:1. 
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(a)  Raw  profile  data 


>  ‘  I  I  I  I  I  I  I  I  I  II  I  I  I  I  I  I  II  I  I  I  I  I  I  I  PROGRAM  A/D  I  |  |  |  |  |  |  |  |  1 1 1 1 1 1  |  1 1  |  |  I  I  |  |  |  t  1  M  I  I  i  t  t  I 
j  10  INTEGER  U<l:30373>,V<l;3O373>,U<l:30373>,D,E,F,Fi1* 

20  DIM  Flit  n*»*«<3> 

21  I  <<<<<<<<<<<<<<<<<<<<<  PROMPTS  FOR  OPERATOR  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 


30 

PRINTER  IS  16 

40 

PRINT  PAGE 

30 

PRINT  LIN<2) 

60 

PRINT  "  To  b*gin 

digitizing, 

pr*ss  CONT" 

70 

PRINT  LIN< 1 > 

80 

PRINT  "  To  stop 

digit izing, 

pr*ss  PAUSE" 

90 

PRINT  " 

not*  lin*  nunb«r  XXX 

100 

PRINT  * 

tvp*  in:  XXX  GOTO  i 

110 

PRINT  “ 

press  STORE" 

120 

PRINT  " 

press  CONT" 

130 

PAUSE 

131  •  <<<<<<<<<<<<<<<<<<<  INITIALIZE  A/D  CONVERTER  >>>>>>>>>>>>>>>>>>>>>>>>>>>>> 

140  OVERLAP 

130  OUTPUT  7, 1; "SN; AC,  1 !  " 

160  OUTPUT  7,l;"MB,O,82;RP,0;RP,80;WB,8;AI,l,l,l;NX;CB;NX!" 

161  !  <<<<<<<<<<<<<<<<  READ  THE  SAW  PROFILE  DATA  FROM  THE  A/B  >>>>>>>>>>>>>>>>>> 

170  FOR  0-1  TO  30375 

180  ENTER  7,1  BFHS  2  NOFORMAT ; <J<0> 

190  NEXT  D 

200  FOR  E-l  TO  30373 

210  ENTER  7,1  BFHS  2  NOFORMAT: V<£> 

220  NEXT  E 

230  FOR  F-l  TO  30373 

240  ENTER  7,1  BFHS  2  NOFORMAT ; W<F > 

230  NEXT  F 

251  !  <<<<<<<<<<<<<<<<<<<  STORE  RAH  PROFILE  DATA  ON  DISK  >>>>>>>>>>>>>>>>>>>>>>> 

260  IF  E-0  THEN  Number  of  f i l*s-l 

270  IF  <E< >0>  AND  <F»0)  THEN  Number  of  f i Us-2 

2S0  IF  F<  >0  THEN  Nu»b*r  of  files-3  ~ 

290  D-D-l 

300  E-E-l 

310  F-F-l 

320  PRINT  PAGE 

330  PRINT  “ NAME " ; Nunb*r_of  f  i Its; “FILES. " 

340  IF  Nunb*r  of  fil*s-l  THEN  INPUT  "  Fil*  nan*  ",Fil*  na»*S<l> 

330  IF  Nunb«r~of~f i 1*s-2  THE  I  INPUT  2  Fil*  t.  *n*s",  FiT«_nair.**<  1  >,  Fi  l*_n 

an«f(2)  ” 

360  IF  Nunb*r_of_f i l«s-3  THEN  INPUT  "  3  Fil*  nuts' , Fi t*_n«b*8< 1 ) , F i I *_n 

«n*f<2>,Fi1*  n«n«8<3) 

370  PRINTER  IS  0 

380  CREATE  Fil*  n*n*8<l>,0,4 

390  ASSIGN  81  TO  Fil*  n*n«*<l> 

400  REDIM  U<l:D> 

410  PRINT  ai;u<*> 

420  PRINT  "  ";D; "POINTS  IN  FILE  "S,Fil*  .)*»•«<  1) 

430  IF  Number  of  fiUs-1  THEN  358 
440  CREATE  FiT*  nan«*<2>,E,4 

430  ASSIGN  82  TO  Fil*  n*n*«<2' 

460  REDIM  V<l:E> 

470  PRINT  82;V<*> 

488  PRINT  *  “  J  E;  "POINTS  IN  FILE  “tFil*  noitS(2) 

490  IF  Nunb«r  of  fil***2  THEN  330 
380  CREATE  FiT*  n*n*t<3>,F,4 

310  ASSIGN  83  TO  Fil*  n«n«*<?> 

320  REOIH  H<l:F> 

530  PRINT  83;H<»> 

540  PRINT  "  ";F;"P0IMT3  IN  FILE  "tFil*  nsn**^* 

330  END 
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le 

26 

21  I 

30 

40 

30 

60 

70 

S* 

60 

90 
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llllllllllllllllllllll  "OGRAH  -Plot-  lllllllllllllllllllllllllllllllllil 
INTEGER  Z< l ! 32767), PI etur«<0! 16300> ,B,C,D,E,M,Ditt,8,Axts 
<<<<<<<<<<<<<<<<<<<<<  PR0HPT8  FOR  OPERATOR  >>>>>>>»»»»>»»>»»>>>>>>> 

PRINTER  IS  16 
PRINT  PAGE 

PRINT  •  This  program! * 

COPIES  A  DATA  FILE  FROM  DISC  INTO  MEMORY* 
INVERTS  THE  RECORD,  AND  REMOVES  THE  2ER0E 


PRINT 

PRINT 


*LCHR6<255:>1<* 

“iCHR*<235>t* 


PRINT 

PRINT 


*LCHR6<2S3)t* 

“LCHRS<235>t" 


PURGES  THE  OLD  RECORD  FROM  DISC* 

CREATES  A  BOAT  FILE  UITH  THE  SANE  NAME  ON 


THE  SANE  DISC* 


100 

PRINT  " 

110 

PRINT  “ 

120 

PRINT  * 

130 

PRINT  * 

140 

PRINT  LIN<1>;* 

130 

PRINT  " 

160 

PRIHT  LIN<1)J* 

ly,* 

170 

PRINT  “ 

34>fc 

*X“tCHR*<34>t"  C 

100 

PRINT  “ 

190 

L INPUT  "  Profi 

200 

LINPUT  *  Graph 

210 

PRINT  PAGE 

220 

PRINT  LIN<10); 

*LCHR»<235>1*  COPIES  THE  NEW  FILE  ONTO  A  DISC* 
*ICHR6<236)S>CHR$(  236  )LCHR6<236) 

"LCHR*<2!3>8."  CREATES  A  PLOT  OF  THE  HEM  FILE  “ 
*LCHRS<233>L"  STORES  THE  PLOT* 

1.  Nwi  the  record  <  £  characters  :  device  >,  * 
press  CONT . * 

2.  Name  the  file  to  receive  the  PLOT,  <  FOR  A  PLOT  ON 


CALL  THE  PLOT  ”ICHR*< 


)* 


press  CONT-. 
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DATA  riOU  BEING  READ  FROM  "tProfile  nue« 

<<<<<<<<<<<<<<<<<<<  READ  RAM  PROFILE  DATA  FROM  DISK  >>>>>>>>>>>>>>>>>>>>>> 


230 

ASSIGN  «1  TO  Profile  nattl 

240 

ON  END  61  GOTO  280 

230 

FOR  Dim-1  TO  32767 

260 

READ  61 S Z<Di a) 

270 

NEXT  Di« 

280 

IF  Dia>-344  THEN  330 

290 

PRINT  PAGE 

300 

PRINT  LIN< 10) 

310 

m 

PRINT  "Record  is  too  small  <*: 

Dim: "points  >  to  use  locating 

320 

STOF 

330 

PUNT  PAGE 

340 

PF'MT  “  FILE  "LProfilt  r.a,n«*S." 

CONTAINS"; Dim; “POINTS* 

330 

PRINT  LIN< 10) J "  FILE  NON  BEING 

INVERTED* 

360 

MAT  Z-Z#<-1>  !  INVERT  THE 

RECORD  HERE 

361 

370 

380 

390 

400 

410 

420 

430 

440 

430 

460 

470 


<<<<<<<<<<<<  REMOVE  STATUS  8ITS  FROM  RAH  PROFILE  DATA  >>>>>>>>>>>>>>>>>>>> 

C-0 

PRIHT  PAGE 

PRINT  LIN<10),*"  ZEROES  NOW  BEING  REMOVED" 

FOR  C-CM  TO  100 
IF  Z<C>»0  THEN  430 
NEXT  C 

IF  <Z<C»81)«0)  AND  <Z<CM62>«0>  AND  '  G'.C*243>«0>  THEN  480 
IF  CO101  THEN  470 

PRINT  “ZERO  LOCATING  ALGORITHM  HAS  FAILED* 

STOP 

GOTO  400 


480  D*0 

490  !  C  IS  THE  LOCATION  OF  THE  FIRST  ZERO 

300  FOR  E>1  TO  Die 

310  IF  <E«C>  OR  <<£-C>  MOD  S:=0/  THEN  340 

320  D«DM  „ 

330  Z  D  •  “Z<  E  > 

340  NEXT  E 

330  N-D  1  N  IS  THE  NUMBER  OF  POINTS  REMAINING 

360  REDIM  Z<l:M» 

361  i  <-.'.-,<<<s<<\<<<<<  PE-STOPE  FAN  RSTFluE  DATA  ON  DISK.  >>*:.•>>'>>>>>,  >'*>>>>  T ' '  > 

370  CALL  R«di»_and_itor«- N,Z- * ',Fr  jf  l 
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4  v„ 


« nw  m mgummfrmtg&M 


371 

380 

390 

600 

610 

620 

630 

640 

630 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

830 

851 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 


I  <<<<<<<<<<<<<<<<<<<<  CREATE  THE  RAW  DATA  GRAPH  >>>>>>>>>>>>>>>>>>>>>>>>>>> 
PLOTTER  IS  13,  “GRAPHICS'' 

IF  Gr4ph  nuiM'X"  THEN  610 

FCREATE  Graph  n *«.*»,  INT<  16381»2.'236>*INT<  16381*2x63536>*3 
GRAPHICS 
GCLEAR 

FOR  B-l  TO  N 

XOV12. 5-120*INT  <  *'B-1  >'  1300) 

Xaxi s*93-4. 75* INT  <B-1)/1500> 

Y“Xaxi s*Z<B>x600 
PLOT  X, Y 

IF  X«120  THEN  PLCT  X,Y,w 
IF  X-120  THEN  MOVE  0,Xaxis 
NEXT  B 


NUMBER  THE  ROHS 


FOR  Axis»l  TO  20 

Xaxi s«95. 4-4. 75*<rtxi s-1 > 

MOVE  121 , Xax i s 
PEN  1 

CSI2E  15x4.54 
FOR  Q«1  TO  20 

MOVE  120+Qx3,X*x!i 
LABEL  “|" 

NEXT  Q 

PEN  -1 

HOVE  120, Xaxi s 
CSIZE  2.2, .5 
LABEL  Axis 
NEXT  Axis 

<<<<<<<<<  STORE  THE  GRAPH  OH  DISK  AND  PROVIDE  ONE  HARD  COPY  >>>>>>>>>>>>>> 

IF  Graph_na**$““X"  THEN  930 

Pictur*<0)«-1 

GSTORE  Pictur*<*> 

FPRINT  Gr  aph_nai»«$,  P 1  cturt<*:> 

PRINTER  IS  0 

PRINT  LIN<5>;“  This  qraph  of  "tChf *<34 >&Pt of i 1 *_n*m*F&CHRF< 34 > 

PRINT  "  has  b«en  stereo  in  &CHR*<34 > LGraph  name *SXHR* <34 ) 

PRINT  LIN<  1 ) ;  “  "&CHR*<34  V:  Prof  i  )  «_na«ne#&CHR*< 34 >8."  cont ai  ns" ; N;  "poi  nt  s.  “ 

DUMP  GRAPHICS 
PRINT  LIN<5> 

END 

SUB  Redim  and  stor *< INTEGER  fi,  INTEGER  Y<*>,F*> 

'NTEGER  X~1:m7 

MAT  X-Y 
PURGE  F* 

FCREATE  F*,  INT<M*2^256»*INT<,1*2x65536)*3 
FPRINT  F*,X<*> 

SUBEND 
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'  iiimtmmimmiiii  "<****  immmmiitimimmimim 

26  • 

30  INTEGER  Za:32?67>,2window<ll61>,Z*dit <t:6t>,«,*,C,D,E,F,6,H, l,Nz,NdSCon,D 
•It*, In  record, XI, X2,  Begin  or  end.X  o,X  Mid* 

40  CON  INTEGER  P i c <8: 16380~, Xval u# , Yval u* 

30  REAL  X, Y 

31  !  <<<<<<<<<<<<<<<<<<<<<  PROMPTS  FOR  OPERATOR  >»»»»»»»»»»>>»»»»> 

60  PRINTER  IS  16 

70  PRINT  PACE 

80  PRINT  "  This  PROGRAM  permits  editing  of  edit  windows  in  *  fil*  contain 

ing  “ 

90  PRINT  "  discontinuities." 

100  PRINT  LIN(l) 

110  PRINT  "  Resulting  corrected  windows  ere  stored  in  e  tep*  file  for" 

120  PRINT  “  us*  b,  PROGRAM  Dscon* 

130  PRINT  LIN< 1 > 

140  PRINT  "  ENTER  THE  FOLLOWING  INFORMATIONS" 

130  PRINT  "  “S.CHR*<223>»."  name  of  the  file  conteining  the  PROFILE" 

160  PRINT  "  " LCHR*  <228  >*. "  name  of  the  fil*  conteining  the  GRAPH  of 

the  profile" 

170  PRINT  LIN(1> 

1S0  PRINT  "  PLACE  A  TAPE  CARTRIDGE  IN  :T14,  end  name  the  fit*  to  record  ed 
it  HINDOUS." 

190  Profit***"" 

200  Graph**"" 

210  Edit**"" 

220  LINPUT  "  PROFILE  fil*  name  .  press  C0NT",Prof i let 

2°0  IF  Profit**-""  THEN  220 

2-.-0  LINPUT  *  GRAPH  fil*  name . press  CONT", Graph* 

230  IF  Graph**""  THEN  240 

260  LINPUT  "  WINDOWS  fit*  name  .  press  CONT", Edit* 

270  IF  Edit*C7, 101<>“ST14"  THEN  2t,0 

230  PRINT  PACE 

290  PRINT  LIN<8> 

308  PRINT  Edit**."  now  being  created.  Takes  about  43  seconds." 

301  •  <<<<<<<<<<<<<<<<<  CREATE  TEMPORARY  CORRECTIONS  FILE  >>>>>>>>>>>>>>>>>>>>>> 
310  CREATE  Edit*, 6101, 4  !  HOLDS  61  INTEGERS  PLUS  X  o  FOR 

320  ASSIGN  »1  TO  Edit*  !  EACH  OF  UP  TO  108  CORRECTIONS  C <61*180)*4 1*4) ) 

330  PRINT  PAGE 

340  PRINT  LIN<8) 

350  PRINT  "  READY  TO  READ  “LProf i te* 

360  PRINT  L IN< 1 ) 

370  PRINT  "  INSERT  CORRECT  DISC .  then  press  CONT." 

380  PAUSE 

390  FREAD  Prof i 1 *«, 2<*> 

400  Nz«R0U<2) 

410  PRINT  PACE 

420  PRINT  LIN<8) 

430  PRINT  «  READY  TO  READ  "IkCraph* 

440  PRINT  LIN< 1 ) 

430  PRINT  "  INSERT  CORPECT  CISC .  then  press  CONT. « 

460  PAUSE 

470  FREAD  Greph*,Pic<* * 

480  PLOTTER  IS  13, “GRAPHICS" 

490  !  BEGIN  TREATING  NEXT  DISCONTINUITY 

300  Next  d:  l 

510  X  as i de*Xl 

320  CALL  Locate  dsconO't  ,X2>  !  FIND  ThE  DISCONTINUITY 

330  IF  <X  asTde-0)  AND  <M  *1>  THEM  660 

546  IF  X1>X  aside  THEN  660 

350  PRINT  PAGE 

560  EXIT  GRAPHICS 

370  PRINT  LIN'S/ 

580  PPINT  "CUPSO*  SHOULD  BE  MOVED  FORWARD  TO  CONTINUE  EDITING." 

390  BEEP 

60O  PRINT  L1H> 1  • 
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610  PRINT  -  to  CORRECT 

0  pro**  CONT. * 

628  PRINT  *  to  RCCEPT 

0  pro**  CONT,* 

630  INPUT  «*,OqP 

<40  IP  Qq«**-*  THEN  920 

<30  IP  Qq«>**‘  THEN  660 

«0  IP  XI  >•!  TNEN  <90 

670  Xt-t 

<00  X2-61 

<90  IF  X2<«ROU<Z>  TNEN  720 

700  X2«R0N<Z> 

710  XI-X2-60 

720  FOR  B«X1  TO  X2 

730  IP  B«X1  THEN  I»0 

740  !•!♦! 

790  Zwindow< I 

760  NEXT  B 

770  X  o-X! 


to  CORRECT  tho  location,  ontor  'tCHR9(34>t*-*tCMR*(34>t*  on 


to  RCCEPT  tho  toeotion,  ontor  ■0CHR8<34>l*o*iCHR»O4>fc*  on 


IE  CERTRIN  TO  STRY  UITNIN  SUBSCRIPTS  ! 


!  Zuindovi'. •>  i  *  not  ro-writton  for  o  window. 
I  SETS  RSI HE  ONLY  ONE  X  TO  LOCRTE  THE  UINDOU 


780  Noxt  try:  HRT  Zodi t *Zwi ndow 

790  ! 

791  !  <<<<<<<<<<<<<<<<<  BRAN  LARGE  GRAPH  OF  EBIT  WINBOW  >>>>>>>>>>>>>>>>>>>>>>>> 

800  Ronow:  GRAPHICS 

818  GCLERR 

820  LINE  TYPE  3 

830  HOVE  0,80 

840  PLOT  0,80 

890  PLOT  120,80 

860  HOVE  0,20 

870  PLOT  0,20 

880  PLOT  120,20 

890  LINE  TYPE  1 

900  HOVE  0,30 

910  PLOT  0,90 

920  PLOT  120,30 

930  FOR  C»1  TO  61 

940  X*<C-1 >*2 

930  Y-50*Z«dit<CV'50 

960  IF  C»1  THEN  HOVE  >.,Y 

970  PLOT  X, Y 

980  NEXT  C 

990  Chooso  typo:  PAUSE 

1000  ~  EXIT  GRAPHICS 

1001  !  <<<<<<<<<<<<<<<<<  CHOOSE  AN  EDITING  PROCEDURE  >>>>>>>>>>>>»>>>>>>>>>>>>> 

1020  PRINTER  IS  16 
1030  PRINT  PAGE 

1040  PRINT  *  BRSED  ON  THE  PLOT,  CHOOSE  ONE  OF  THE  FOLLOHING: " 

1090  PRINT  LIN< 1 ) 

1060  PRINT  "  *tCHR*<240)t •  to  viow  tho  plot  »g»ir>,  ontor  ?* 

1070  PRINT  *  "tCHR I <240 to  »kip  thit  locotion,  ontor  /• 

1080  PRINT  LIN<1> 

1090  PRINT  *  DISCONTINUITY  TYPES  :* 

1100  PRINT  LIN< 1 > 

1110  PRINT  "  "fcCHRf <240. t "  BEGINNING-  OR  END-OF-RECORD,  ontor  R* 

1120  PRINT  *  “tCHR#<240  .  STRAIGHT-SLOPE  REHOVAL,  ontor  S“ 

1130  PRINT  «  "tCHRP<240''4 *  POINT-BY-POINT  FIX,  ontor  P* 


1130  PRINT  * 
1140  PRINT  * 
1130  Typo*-"’ 
1160  INPUT  " 


*t  CHRP* 240>»"  TIE  <  NITN  STRAIGHT  LINE  )  ,  ontor  T“ 


?  or  r  or  R  or  S  or  P  or  T 


thon  C0NT*,TypoP 


1170  IF  <TypoP<  >*?*  ’  AND  <Tyt«f<  Al,0  TypoP 

poPO"P")  AND  < TyptPv  >  *T"  >  THEN  1160 
1130  IF  TypoP»"R*  THEN  CALL  Ptcord*.!  o.  ’odi*<*>* 
1190  IF  TypoP»*S*  THEN  CALL  S '  cf.«  I ,  Cou  •«*■»•  > 
1200  IF  TypoP«"P"  THEN  CALL  P: in* * lodi t t * . ■ 

1210  IF  TypoP*"T"  THEN  CALL  S **p. <0.r*dit < * > > 
1220  IF  TypoP<  THEN  1250 
1230  GRAPHICS 

1240  GOTO  Chooto  •  vp« 

1250  IF  TypoP"*'*  THEN  No<t  o“ 

1260  PRINT  PAGE 


Al,8  TypoPO'R")  AND  <TypoP< >“S“ >  AND  <Ty 
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127*  MINT  LIN<9)  \ 

!2M  MINT  •  IF  SATISFIES  H1TH  EDITING  ON  THIS  NINBOH,  ent ©r  \ 

I2M  -Ft  I  NT  TO  M-CB|T<4«-4N9|  MtlN  FOt  TtMt  NINBON,  enter  -•  \ 

13M  MINT  '  TOT  tONTlNNfc/TNIi  HINbOM  UITN  ANOTHEA  TVFt  OF  COMICTION,  enter 

I  I2IO  iMitrutlenM*' 

1320  INFUT  ■  ♦  or  -  or  s', Sot tsf set  ion* 

1330  IF  <9«t 1 sf act  t  on*< >“♦“ )  ANB  <Set i *f set i on«< >■-■ >  AN8  <Soti»foction»<>*>'*> 
i  THIN  1320 

134#  IF  SattsfactionOe"**  THIN  Set  aside 
1330  IF  Satisfaction*-*-*  THIN  Next  trv 
134#  IF  Sat i  sf  ac  t  i on*- * / *  THIN  Rene* 

1341  I  <<<<<<<<<<<<<<<  HAITI  IBIT  HI  NOON  ONTO  TAFI  CAATAIBCE  >>>>>>>»>»»»»> 
1370  ltt_uiO!  I  Set  uidt  <  12  points  >. 

13AA  FAINT  »i;X  o 

1390  FAINT  Ol;Zedit<#> 

1400  GOTO  Next  d 

I4At  i  <<<<<<<<<<  SUBROUTINE  FOF  RICOABINC  ICGINNINC  OA  INB  OF  FILO>>>>>>>>>»> 
1410  BUB  Record< INTECIA  Xzero, INTEGER  2<e>> 

142#  GRAPHICS 

143#  INTECIA  X1,N2,N  1,2  wtndow( 1 :41 > , A 
144#  AIM.  U,V 

148#  CSI2C  I9'4.94 

144#  NOV I  2#,9# 

147#  LABEL  ‘RSCOAB  M  INB  M  START* 

14##  NOVI  20, #9 

149#  LABEL  *<  point  ,  then  tJ  >* 

If##  M1NTIA  1#,99,# 

19#t  i  <<<<<<<<<<<<<<  M0V1BI  CUASOA  TO  LOCATE  BOF  or  EOF  >>>>>>>>>>>>>>>>>>>>>> 

191#  BICITIZI  XI, VI 
t#2#  BICITI2C  X2.V2 
193?  IF  Xl<2  THEN  Xl*2 
1**  IF  XI >11#  THEN  XI* 122 
I99C  X  l*X«ere*IMT«Xl/2>-l 

IBM  TF  < V2>#9>  AN#  <X2>94>  AND  <X2<57>  THEN  Nsrk-1 
t#7#  IF  <V2>#9)  ANB  <X2>7C>  ANB  <X2<79>  THEN  Hark-2 
l«M  IF  ltark-1  THEN  2  wi  ndow<  INT<X1^2>  >—32747 
199#  IF  Nark -2  THEN  2~windew<INT<Xl< 2>>-32747 
I###  C812S  2.9 
1«1#  NOVS  Xl-I#,#2 

I #2#  IF  Nark-1  THEN  LABEL  *  ENB  I  ‘LVALFCX  l> 

1C3#  IF  Nork -2  THEN  LABEL  "START  |  *lVAL»<x”l > 

144#  C#'TE  13  4.34 
149#  LINE  TVM  3 
144#  HOVE  XI, #1 
147#  FLOT  XI, #1 
14##  FLOT  XI, 19 
149#  LINS  TYFE  1 
179#  PAUSE 
171#  MIHTEA  IS  14 
172#  MINT  FACS 

173#  FAINT  L1N<9>;*  If  tstisfiod  with  tho  correction,  enter  ♦" 

174#  FAINT  LIH<1>;*  To  re-do  the  correction,  enter 

173#  #©••*“ 

174#  INFUT  *♦  ,  -*,Bo» 

177#  IF  <8oS<  >*♦*)  ANB  <B©*<>-  )  THEN  1740 

17##  IF  B:t-"-*  THEN  1610 

179#  2<INT<X1'2>)>2  window*  IMTc:i-'2)  )  '  CHANGE  ONLY  THE  VALUE  OF  ONE  Z<*> 

1###  SUBEXIT 

t##l  i  <<<<<<<<<<<<<<<<<  GAAFH  OF  THE  COPRECTEB  EBIT  HINBOH  >>>>>>>>>>>>>>>>>>>> 

l#l#  CCLEAR 
1620  LINE  TYPE  3 

1930  HOVE  9,60  I 

1840  FLOT  0,69  , 

1990  FLOT  120,90 

1040  HOVE  0,20  I 

11670  FLOT  0,20  { 

11860  FLOT  120,20  i 

|1890  HOVE  0,50 

1900  LINE  TYFE  1 

1910  PLOT  0,50 
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192*  PLOT  129f99 

1939  FOP  9-1  TO  61 

1949  U-<9-l>«2 

1939  V-S9*2<9V90 

1999  IF  9-1  THEN  HOVE  U,V 

1979  PLOT  U,V  ■> 

1999  NEXT  9 
1999  GOTO  1469 
2999  SUBEND 

2991  •  <<<<<<<<<  SUBROUTINE  FOP  SLOPE  OR  STRA1 GHT-LINE  CORRECTION  >>>>>>>>>>>>>> 


2910  SUB  S I op«< INTEGER  Slop*  or  tit, INTEGER  Z'#>> 

2029  INTEGER  XI , X2, X3, Y1 , Y2, Z_w7ndow< 1 : 61 ) , A 

2930  REAL  U,V, Slopt  fix  1, Slope  fix  0,Sttp  walut  l,Sttp  oalut  0,Y3,Rtal  z 
2940  GRAPHICS 

2030  CSIZE  13  '4. 34 

2060  HAT  Z  wi ndow*Z  <  VALUES  OF  NEED  NOT  BE  CHANGED  OR  NOTED 

2070  HOVE  10,90 

2080  IF  Slop#  or  tit-1  THEN  LABEL  "STRAIGHT-LINE  SLOPE" 

2090  IF  SI opt~or~t i t-0  THEN  LABEL  “STRAIGHT-LINE  TIE  “ 

2100  POINTER  To, 35, 8 
2110  DIGITIZE  XI, Y1 
2120  DIGITIZE  X2.Y2 

2139  DIGITIZE  X3.Y3 

2140  IF  X3-X2  THEN  2220 

2150  PRINT  PACE 

2160  EXIT  GRAPHICS 

2170  PRINT  LIN<8; 

2180  PRINT  ■  LAST  2  X'*  MUST  BE  EQUAL." 

2190  PRINT  “  Start  ovtr  by  pressing  CONT. “ 

2200  PAUSE 

2210  GOTO  2100 

2220  Xa-INT<Xlx2)*l 
2230  Xb-INT<X2/2)*1 


22*0  Ptal  z-30*<Y3-5O) 

2250  IF  Xb>Xa  THEN  2330 
2260  EXIT  GRAPHICS 
2270  PRINT  PAGE 

2280  PRINT  LINC5) ; "  Sine*  X<a>*";Xa;"  and  X<b >»" ; Xbi " ,  tht  slopt  has  tht  wrong  s 
i  gn.  “ 

2290  PRINT  "  PLEASE  REPEAT,  and  make  X<b>  >  X<a>." 

2300  BEEP 
2310  PAUSE 
2320  GOTO  210O 

2330  Sttp_val ut_l*R*a1_z-Z_wi ndou(Xb) 

2340  Sttp_va1  ut_0«Rtal_z-Z_wi  ndo'j<Xa> 

2350  Slop*_f  i  x_l*St*p_val  u*_lx<;:b-X*> 

2360  Slopt  fix  0-Sttp_valu«  0x<Xb-Xa> 

2370  FOR  A-l  TO  61 
2380  IF  A< -Xa  THEN  2440 

2330  IF  <S1opt_or_ti#-l)  AND  (A>Xa>  AND  <A<Xb>  THEN  Z_wi ndew< A>-Z_wi ndo w< A)*S1 op 
t_f  jx_l#(.A-Xa> 

2400  IF  <S1 op*_or_t i e-0>  AND  <  A  >'.;*>  AND  <A<Xb)  THEN  Z_wi  ndow<  A)-Z_wi  ndowtXa' ♦  $  1  o 
pt  fix  0*<A-Xa>  ~  ” 

24?0  IF  < A>Xa>  AND  <A<Xb>  THEN  LINE  TYPE  3  '  DOTTED  LIME  IN  FIX  AREA 
2420  IF  Ax-Xb  THEN  Z_wi ndow< A >*Z  wi ndow< A) +St *p  value  1 
2430  IF  A>-Xb  THEN  LINE  TYPE  1 
2440  U-< A-l >*2 


2450  V-Z  uindow<A-x3O*50 

2460  IF  A-l  THEN  HOVE  U,V 

2470  PLOT  U,V 

2480  NEXT  A 

2490  PAUSE 

2300  PRINT  PAGE 

2510  PRINT  L I N < 5 


2320 

PRINT  " 

TO 

ACCEPT 

THE  CQf SECTION,  tn'tr  ♦ 

a  a  .prtll 

CONT" 

25  30 

PRINT  " 

TO 

CHANGE 

THE  COPPf CTIQN,  «n»tr  - 

a  a  .prtS* 

CONT" 

2540 

If#*"" 

2350 

INPUT  " 

♦ 

or 

-".If* 

2360 

IF  (If*. 

>  ♦ 1 

")  AND 

•  If#  "-">  THEN  2530 

2370  IF  If#-"-"  THEN  2600 


28 


2989 
2998 
2991 

2898 

2(19 

2(29 

2(39 

2(49 

2(99 

2((9 

2(79 

2(89 

2(99 

2799 

2719 

2729 

2739 

2749 

2799 

27(9 

2778 

2788 

2799 

2899 
2891 
2818 
2828 
2838 
2848 
2898 
28(9 
2878 
2889 
2t90 
2908 
2918 
2929 
2938 
2949 

2990 
29(8 
2978 
2989 
2998 
3000 
3818 
3829 
3838 
3849 
3999 
39(0 
3870 
3980 
3098 
3180 
3110 
3120 
3130 
3140 
3190 

Z 

31(9 

3170 

3180 

3190 

3200 

3210 

3220 

3239 

3240 


HOT  Z-Z  window 

SUBEXIT- 

1  <<<<<<<<<<<<<<<<<  GRAPH  OF  TMC  C0RRECTC8  t»IT  MIHBOW  »>>>>»»»»»>>>> 

GRAPHICS 

GCLCOft 

LINS  TYPS  3 

HOVS  9,89 

PLOT  9,99 

PLOT  129,80 

HOVE  9,20 

PLOT  9,20 

PLOT  129,29 

LINS  TYPE  1 

HOVE  9,90 

PLOT  9,99 

PLOT  129,99 

FOR  AM  TO  (1 

U-(A-1 >*2 

V-90*Z<A)'99 

IF  A-l  THEN  HOVE  U,V 

PLOT  U,V 

NEXT  0 

GOTO  29(8 

SUBEND 

1  <<<<<<<<<<<<  SUBROUTINE  FOR  POINT-BY-POINT  CORRECTION  >>»»>>>>>>>>>>>>> 

SUB  Point (INTEGER  2<*>> 

INTEGER  XI , X2, Z  u i ndow< 1 : (1 > , Z  plot ( 1 :(1 >,0,C,B 

REAL  Y1,Y2,U,V,Z  r*al,D#lt* 

CSIZE  IS'4.94 

Start:  GRAPHICS  I 

HAT  2  window-Z 

Noxt :  ”  1 

HOVE  18,90 
LABEL  "POINT  FIX- 
HOVE  2,9 

LABEL  "<  up  or  down.” 

LINE  TYPE  3 
HOVE  2,100 
PLOT  2,108 
PLOT  2,0 
LINE  TYPE  1 

Horc:  I 

POINTER  10,99,8 
DIGITIZE  XI, Y1 
DIGITIZE  X2.Y2 
IF  X2>-2  THEN  3098 
X1-X2-0 

Dol t»"<Y2-Yl )*98  I  PROVIDE  FOR  HOVING  ENTIRE  MINDOH  UP  OR  DONN 
GOTO  Plot 

IF  X2-X1  THEN  None* 

EXIT  GRAPHICS 
PRINTER  IS  1( 

PRINT  PACE 

PRINT  LIN<9>,"  X'l  mu**  bt  «qu*t . . .  prt»*  CONT" 

PAUSE 

GRAPHICS 

GOTO  Nor*  •  TRY  DIGITIZING  AGAIN 

Hone*:  • 

Z  ro«l-S8*<Y2-99> 

D#lt*-Z_r**1-Z_window(<X2*2W>  I  THE  CORRECTION  FOR  THIS 

Plot:  i 

GRAPHICS 

HAT  Z  plot -Z  window  •  SET  A  TEHPOPARY  PLOTTING  VECTOR 

FOR  C-l  TO  (1 
U-<C-t>*2 

IF  U'-XI  THEN  LINE  TYPE  3  »  DOTTED  LINE  AFTER  CORRECTION. 

IF  U>-Xl  THEN  Z  plottO-Z  plot <C)*D*lta 
V-Z  plot  (O-  50*?0 
IF  C-l  THEN  HOVE  U,V 


29 


* 


3250 

3260 

3278 

3268 

3298 

3368 

3?16 

3328 

3338 

3348 

3398 

3368 

3378 

3388 

3390 

3488 

3418 

3428 

3438 

3440 

3458 

3451 

3460 

3470 

3488 

3498 

3500 

3510 

3528 

3538 

3548 

3558 

3568 

3578 

3588 

3590 

3688 

3610 

3620 

3638 

3648 

3658 

3660 

3678 

3688 

3698 

3708 

3718 

3728 

3738 

3731 

3748 

3750 

3760 

3778 

3788 

3790 

3800 

3810 

3920 

3830 

3840 

3656 

3860 


PLOT  U,V 
NEXT  C 

LINE  TYPE 

PAUSE 

EXIT  GRAPHICS 
PRINTER  IS  16 
PRINT  PAGE 

PRINT  LIN< 5) ; “ 

PRINT  LIN<2>;“ 

PRINT  LIN<2>;“ 

PRINT  LIN<3>;“ 

Sat**-" 

INPUT  "“.Sat* 

IF  <Sat*<>"+“>  AND  <Sat*<>“-“>  AND  ( Sat  *<>•»“  >  AND  <Sat»< >■✓" >  THEN  3370 


If  this  correction  is  Acceptable,  enter  ♦* 
If  this  correction  is  NOT  Acceptable,  enter  - * 
To  begin  with  the  or i gi nal  window,  enter  s‘ 
If  all  corrections  are  complete,  ei^ter  *“ 


IF  Sat*«“*“ 
IF  Sat*»"-“ 


THEN  Final 
THEN  Minus 


Plus: 


Minus: 


IF  Sat*»“'“  THEN  Minus 
IF  Sat*»“+“  THEN  Plus 

MAT  Z_window*»Z_plot 


1 


THE  PLOT  VECTOR  IS  NOM  THE  EDIT  HINDOU 
CONTINUE  BY  FIRST  RE-DRAHING  THE  PLOT 
•  <<<<<<<<<<<<<<<<<  GRAPH  OF  THE  CORRECTED  EDIT  HINLOH  >>>>>>>>>>>>>>>>>>>> 

GCLEAR 
GRAPHICS 
LINE  TYPE  3 
HOVE  0,80 
PLOT  0,80 
P LOT  120,80 
MOVE  0,20 
PLOT  0,20 
PLOT  120,28 
LIME  TYPE  1 
MOVE  8,50 
PLOT  0,50 
PLOT  120,58 
FOR  AM  TO  61 
UMA-1>*2 

IF  <Sat*a“-">  OR  CSatM-e-*)  THEN  V«50*Z  wi ndow<A)^50 
!  LINE  3248  SHOWS  Z  window  TO  BE  THE  PRESENT  WORKING  EDIT  WINDOW 

IF  Sat#»«/“  THEN  V-58+Z<A)/50  !  RETURN  TO  ORIGINAL  Z  window 

IF  AM  THEN  MOVE  U,V 
PLOT  U.V 
NEXT  A 

IF  Sat *■"/■*  THEN  Start  !  TEGIN  .'GAIN  BE  RE-LOADING  Z<*> 

IF  <Sat«»--")  OR  <Sat**“t" )  THEI.  :«xt 
Final:  1  ACCEPT  LOIT  WINDOW  • 

MAT  Z_window»Z_plot 
MAT  Z«Z  window 
EXIT  GRAPHICS 

SUBEND 

I  <<<<<<  SUBROUTINE  FOR  LOCATING  EDIT  WINDOW  USING  SMALL-SCALE  PLOT  >>>> 

SUB  Locate_dscon< INTEGER  Xzero, INTEGER  XI08> 

INTEGER  Correct ion_at_x 
COM  INTEGER  Pic<*>,XMAlue,  lvalue 
GLOAO  Ptc<*> 

GRAPHICS 

POINTER  Xualue, Y"Alue,2  '  POINTER  IS  PLACED  AT  LAST  CORRECTION  POINT 
!  BEGIN  CORRECTIONS  HERE 
DIGITIZE  Xualue, Yvalue 

Line_nuaber»INT<<l82. 125-Y«-sl ueV4. 75 > 

Correct  ion_at  _x«  1580 *L  i  ne_nuab«r“I499*Xy*lue*-12. 5 
•XieroeCorYec t  i  on_at - 38~ 

Xt08»Correct ion  at  > *38 
SUBEND 


>> 
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>  *  I  II  1  I  II  I  I  I  I  t  I  II  I  I  I  I  I  I  I  I  WOGMH  “ti scon*  |  |  |  |  |||  I  1 1  I  III  I  I  1 1  I  I  I  I  1 1  I  1 1  I  I  I  I  I  !  I  I 

2  <  <<<<<<<<<<<<<<<<<<<<<  PROMPTS  FOR  OPERATOR  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 

PRINTER  IS  16 
PRINT  PAGE;LIN<2) 


10 

20 

30 

40 

SO 

60 

70 

00 

*e 

.00 

110 

120 


PRINT 
PRINT  " 

PRINT  L!N<2> 
PRINT  “ 

PRINT  " 

PRINT  " 

PRINT  LIN<1> 
PRINT  “ 

PRINT  LINO) 
PRINT  " 


This  progra*  uses  th*  edit  windows  recorded  by  program  Edit* 
to  rewrite  *  profit*,  taking  out  di scont inul t Its. " 


ENTER  th*  n't.*  of  th*  f11*“ 

containing  th*  PROFILE 
containing  th*  edit  WINDOWS 


<  *  BOAT  fit*  >," 

<  a  BATA  f 1 1*  >. ■ 


LE  NAME.  >* 


Then  press  CONT." 

<  AS  A  RECORD  IS  COMPLETED,  ITS  NAME  NILL  BE  REQUESTED  AS  A  FI 


130 

140 

ISO 

160 

170 

180 

ISO 

200 

210 

220 


Pf.U«-‘ 
L INPUT  " 


PROFILE  ...  th*n  CONT",P$ 

IF  P*18;l]0"F*  THEN  140 
LINPUT  *  WINDOWS  ...  then  CONT", UR 

INTEGER  Z< 1 : 301O0), Wi ndot< 1 : 106, 0: 61 : , A, Z_si z*,D*1 t  a_z 
INTEGER  No  of  wi ndous , Loop, X, Current  x, Window  x,Next~x,2  Aside 
ASSIGN  II  TO  Ml 
FREAD  P»,Z<*> 

ON  END  HI  GOTO  Continue 
PRINTER  IS  0 

221  !  <<<<<<<<<<<<<<<  READ  EDIT  WINDOWS  FROM  CARTRIDGE  >>>>>>>>>>>>>>>>>>>>>>>>> 

230  FOR  A*1  TO  100 

240  FOR  B*0  TO  61 

258  READ  il ; Wi ndow< A, B)  !  Uindow<a,0>  is  th*  window  st Anting  locAtion. 

260  NEXT  B 

270  PRINT  Window<A,0>; " | 

280  NEXT  A 

281  !  <<<<<<<<<<<<<  BEGIN  CORRECTIONS  (All  automatic)  >>>>>>>>>>>>>>>>>>>>>>>>>> 

>>>>>>>>>>> 


2S0  Continue:  ! 

300  Z_si ze*R0W<2) 

310  No_of__wi  ndows«A-l 

320  Current  xM 

330  !  ~ 


340 

•  AS* 

•  * 

First  window  <  contains  a  DOR  ) 

350 

! 

11.  Hi ndow<  1 , 8)M 

360 

! 

12.  Window! 1 ,0>< ' 1 

370 

!  Cos* 

2: 

Encountering  a  new  window  <  Current  x>Next  x  >  with: 

380 

( 

21.  Window  x*0 

390 

211.  Loop»No  of  windows  (  contains  EOR  > 

400 

• 

22.  Window  xOO  (  overlapping  windows  ) 

410 

« 

221.  Loop-».ic  of  windows  (  contains  EOR  ) 

420 

!  CAS* 

3: 

In  a  window  not  LAST  window<  Current  x<N*xt  x  AND  Window 

430 

1 

Encounter ing  a 

440 

1 

31.  DOR  <  Window (Loop, Window  x«32767  ) 

450 

1 

32.  EOR  <  Window  Loco, Window  **-32767  ) 

460 

•  Cas* 

4: 

Leaving  a  window  <  i.4.  Window  x*61  > 

470 

•  Cas* 

5: 

Betiison  windc.*  i  <  Cur-«n»  -(Next  x  AND  Uindou  x*0  ) 

480 

1 

490  IF  <Curr*nt_xM  •  AND  • windovC 1 , 0>M )  THEN  Cesell 

500  IF  <Curren»_xM >  AND  • Window  1 , 0>  >1>  THEN  C*S*I2 

510  N*'t_point : ! 

5*0  IF  '  Current _x“N*'<t_*  *  HMD  •.  Mirdoi-_  '»0>  THEN  Cas*21 

520  IF  :  Curren*  _x*Ni  AND  '  Wi  ndow_v  >0  >  AND  'Wmdow_  -61)  THEN  Cas«2C 

540  IF  •  Curren*~x- t_<  AND  •  Wmdo"_  ->0.'  AND  <Ni ndow_- <61 1  THEN  Ca**3 

550  IF  Windou_\a£i  THEN'Ca»<4 

560  Ii;  'Current  x  N*  t  •  •  AND  ‘Window  »0 >  THEN  C*s*S 

570  Ckteli: 


31 


580  Loop"! 

599  Next  x"Mind ou(2,9) 

699  FOR  A"1  TO  61 

610  IF  UindouC 1,A)< >32767  THEN  679 

629  Rrint«"“0N* 

630  Mtndou< 1 ,fl)*Hindowv 1 , AM ) 

649  Current_x»A 

658  Mindow_x"A 

669  GOTO  Next  point 

679  NEXT  A 

689  C*sel2l 


!  ASSUME  Next  x  occurs  AFTER  BOR 


698 

700 

718 

728  C*se2l: 


Next_x"Current_x"Mi ndou< 1,8) 
MindowO,  1  >*Mi  ndou<  1 , 2) 

GOTO  Next_point 


7 38  Loop>LoopM 

748  IF  Loop=No_of_ui ndows  THEN  C*se211 

758  Next  x»Wi ndou(Loop+l , 8) 

768  IF  Print««"0FF“  THEN  SkTp2l 

778  Deltx  z*Z<Current  x)*Delt*  z-Wi ndow<Loop, 1 ) 

788  IF  A8S<Mindow<Loop, 1>+Delt*_z >“30080  THEN  878 

790  !  !  ;  i  |  1  M  !  I  !  !  I  I  !  I  M  !  !  !  i  M  M  I  !  M  nTi  !  !  !  M  !  M  !  !  !  I  M  !  !  !  I  !  !  !  !  I  i  I  !  M  !  I 

888  REDIM  Z<1:X>  (_ 

810  CALL  Record_record<Z<*> , X)  !_ 

928  REDIM  2<1:Z  size)  *  ! 

830  IF  MindouCLoop,  1  >-»De1  te_z-.  -30800  THEN  Be  1  t*_z-De  It  *_2  *38000  •“ 

848  IF  Mi ndowCLoop, 1 >*De1 t &  z  30800  THEN  Belts  z«De1t»  z-300O8!_ 

858  X»0  “  !  _ 

868  !  t  m  !  I  !!!!!  I  !l  !  (!!!!!!!!!!  M  !!!!!'!!!  M  !!  M  !!!!  M  !!!!!!!  M  I  !!!  lT 


878 

830 

898  Skip2i: 

988 

910 

920 

939  C*se22: 


X«XM 

Z(X)"Mi ndouCLoop, 1 )*Del«  *_z 
Mindou_x*2 

Current _x*Current_xM 
GOTO  Next_point 


!  DO  NOT  RE-COMPUTE  Deltk  z 


948  Loop»Lcop+l 

958  IF  Loop*No_of_wi ndows  THEN  C*sc221 

968  Next_:<*Window<Loop*l,8> 

970  IF  Pr i nt  S*"0FF"  THEN  Sfcip22 

980  Belt*  z"Hi  r>dow<Loop-l ,  Mi  ndow_x>-»Be  1 1  *  z-Mi ndow<Loop. ’ > 

998  IF  ABSCMi ndouCLoop, 1  >  +  De 1 t *_z><  3O080  TMEN“i088 

1800  !  !  !  H  M  !  II  I  !  !  I  !  M  !  !  !  !  H  M  !  •  H  •  f  Tl  !  f  •  •  !!  !  !  H  !  !  I  !  !  !  i  *  t  I  !  H  •  1  I  !  M  1 

1010  REDIM  Z< 1:X>  !_ 

1929  CALL  Record  record<Z(»),X>  ! 

1939  REDIM  Z<l:Z~size;  C 

1040  IF  Mindow<Loop, l)*Delt*_z<.-3O800  TUEN  Del t *_z«Del t*_ze30003 
1858  IF  Mi ndowtLoop, 1 >+De I t *  z;3006o  THEN  Belt*  z»Delt*  z-38803! 

1060  X-8  ‘  ~ 

1070  !  !!!!!!!!!!!!!!!II!!!!{{!!I!!!!!!!!!!!|H!!M!1!!MI!I!!!!!!!H! 

1080  X«XM 

1090  Z<X)"Mindow>.Loop,  1  )eDelt*  z 

1100  Sk  i  p22l  “  I  DO  NOT  RE-COtlPUTE  D*lt*_z 

1119  Mindow_x"2 

1128  Currer>t_x"Current_xM 

1139  GOTO  Next  point 

1140  C*se21 1 :  “  ' 

1150  D#lt«_z«Z<'Currert_^  >eDel  t  a_z -w i  i«3cv'  Loop,  1) 

1160  GOTO  Go  on  21 

1170  Ckse221!  “  “  I 

1188  De1t»_z"Mindow< Loop-1  ,Umdou_  )»D*lt*_z-Mindou<Loop.  1  > 

1198  Go  on  21 :  “  ~  > 


FOR  AM  TO  61 

IF  Window Loop, A*  -32767  THEN  :Kip__21 
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1220 
1230 
1240 
1290 
1200 
1270 
1200  Skip  211 

I  “  ' 

1290  IF  ABS(Hindow<Loop,A’eDelt a_z><  3OO00  THEN  1380 

1380  I  !  !  !  I  I  !  !  !  !  I  M  !  I  i  M  H  I  !  !  !  I  !  !  !  M  1  •  !  •  >  ■  i  i  1  M  i  i  i  i  I  1  I  i  i  !  M  !  I  !  i  *  *  I  < 

REDIN  2<i:x> 


REDIM  M!X) 

CALL  Record  record<2<*>, X> 
FRINT  PACE  " 

PRINT  LINC5) 

PRINT  “PROGRAM  COMPLETES a 
STOP 


1310 

1320 

1330 

1340 

1350 

1300 


CALL  Record  reeord».2<*>,X> 

RESIN  Z<l:Z_size> 

IF  U i  ndow < Loop,  A >'•’&•  I  t«_z<- 30000  THEN  De1ta_z-Delta_ze-38000 
IF  Mindow(Loop,A)eDelt*  z>30000  THEN  Belt*  z-Del ta~z-300O0 
X-0 


1370  1  M  f  !  j  (  !  I  t  !  1  !  !  •  !  M  M  !  I!  •  M  M  !  M  !  !  !  !  !  !  1  !  !  !  M  !  M  I  !  !  M  !  !  H  I  11  !  !  !  I  I 

1360  X-X+l 

1390  Z<X)-Uindow(Loop,A>eDtlt4  z 

1400  NEXT  A 

1410  RESIN  2-l:X> 

1420  CALL  Recc-rd_record(Z<'0,X; 

1430  PRINT  PACE 

1440  PRINT  LI  IK  5) 

1450  PRINT  “  PROGRAM  COMPLETED*  . 

1400  STOP 

1470  Case3: 


I486  IF  HindouiLoop,Mindou_>0-3276?  THEN  Cais«31 

1490  IF  Mi  ndowt  Loop,  Mi  nd0‘-_x>-- 32767  THEM  Ce»e32 

1500  GOTO  Case  3  continue 


Pr int*-"ON“ 

Mi  ndotKLocp,  Mi  ndcu__x)*Wi  ndouCLoop,  Hi  ndou_xel  > 
De1ta_z-0  ~  ~ 

GOTO  Case  3  continue 


Print«-"OFF" 

Delta  z-8 
REBIN~2(i:X) 

CALL  Record  record<2< *  1 , X' 

X-0 

REDIM  2<l:Z_»ize> 

1030  Ca*e_3_eont inue: 

1040  IF  Print •■"OFF"  THEN  No  print_3 

1050  IF  ABS<Nindow<Loop, Mindow_  ♦D* lta_z><30000  THEN  1740 

1800  I  !!!!!!I!!!!!)!!!!!!!!!!I{!I!*!!I!!!!iT!!!!!!!1!!!!!!I!!!!!!!!!!!!!I 

1670  REDIM  Z(t:X) 

1680  CALL  Record  recor dt2<*>»X> 

1690  PESIM  2ti:2~*ize  • 

1780  IF  Hindow<Loop,Mindow_  >eDeUa_z<-3C’O0O  THEM  Delta_z-Del ta_ze30OO0 
1710  IF  Mindow<Loop,Hindo~  •♦Delta  030000  THEN  Delta~z-Delta~z-3O8O0 

1720  X-0  “ 

1 738  I  » i ! i i i | j  j I | : j i M i i M ! i i i i I M l II i i i i j ii i j i i I j  I  |  It  !  I  l  j  I  |  M  •  !  !  M  1  !  M  !  ! 
1740  x-:*i 

1 750  Z<  V1  -M  t  ndou  < Loop , W i ndow_x ) +De 1 1 a_z 

1760  No_print_3: 

1770  Current  -Current  .el 

1780  MindOu_”-Mi ndow_  M 

1790  GOTO  He  t  point 

1800  Cate4: 

0 

1810  IF  Pr  int  •■■OFF' "  Tnf.  fMr4 

1820  Delta  z-Ht  n«u'  L  .  p.t-l  *Dtlt  i  :*2'Current_x' 

1830  IF  ABitZ-  Current  ..-telt  *  ;  •  v  30vO*’  THEM  1920  “ 

1840  1  i  i  i  i  t  i  i  l  i  i  i  i  M ii i i i T i  i  •  i  i i V i  i  •  i  ■  i  i  •  i  i  i  i  i  i  i  ii i i i i i i i i i i i i i i 


1510  Case3l: 

1520 

1530 

1540 

1550 

1560  Caee32: 

1570 

1580 

1590 

1600 

1610 

1620 
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1850  REDIM  2<1:X>  !_ 

1880  CALL  Record  record',2<*>,X>  !~ 

1870  RCDlil  2<  l:2~S'=*  '  t“ 

1880  IF  2<Current_x>eDel t*_z'-30OO0  THEN  Be  1 1 *_z*De 1 1 *_z*30000 ! ~ 
1880  IF  2<Current_:<>eDel  t*  z>30000  THEN  Delta  z*De1ta  z-30008  •  ” 
1980  ~  X-0"  “  i~ 

1918  !  M  i  i  |  i  M  i  i  i  •  !  M  !  M  M  M  '  !  H  H  !!  !  M  !  H  •  M  !  •  M  !  M  M  M  M  •  !  !  !  •  !  !  T 
1920  X*X*1 

1930  2<X)*Z<Current  z 

1940  Skip4: 

1950  Windou_<*0 

1980  Curr#nt_x*Current_ye  1 

1970  GOTO  Next  point 

1980  Case5: 


1990 

2088 

2010 

2028 

2038 

2040 

2058 

2080 

2070 

2080 

2090 

2108 

2110 

2120 

2138 

2140 

2158 

2151 

2180 

2178 

2180 

2198 

2200 

2210 

2228 

2230 

2248 

2250 

2280 

2278 

2280 

2290 

2300 

2310 

2320 

2330 

2348 

2350 

2380 


FOR  ft*Curr*nt _x  TO  Next 1 
IF  Pr  i  fit  S*“OFF“  THEN  2120 
IF  2 < ft) +Lel  t  30800  THEN  2108 

!  >  !  !  !  !  !  !  !  )  !  !  !  |  !  !  !  !  !  i  |  !  |  !  i  !  i  |  hi  |  |  |  !  {  ;  i  i  ;  ;  |  |  |  |  |  {  i  ;  ;  j  |  i  {  {  |  m  |  ;  | 

REDIM  Z < 1 1 K >  !_ 

CALL  Record  record<2<*> , X)  ! 

REDIM  2<l:Z'size> 

IF  2< ft ) +De  1  t*_z< -38060  THEN  De  I  *  a_z=De 1  t  a_z*30.O00  ! _ 

IF  2<ft)*Del t %_z >30008  THEN  lei t a_z=De 1 t a_z-30G00  ! 

X*0  ' 

!  !  !  !!!!!!!!  H  !  1  !  !  !  !  !  !  !  !  *  1  !  !!!!!!  H  M  !  1  !  1  !!!•!!  H  !  H  !  1  M  !  M  !  H 

);*:<♦  l 

2<X)*2(fl  > ♦De 1 «  a  Z 

NEXT  ft 

Currerit_x=-Next_;<. 

GOTO  Next  point 

END 

!  <<<<<  SUBROUTINE  FOR  STORING  THE  DftTft  UtTH  DISCONTINUITIES  REMOVED >> 
SUB  Record  record*!  INTEGER  2  inf *), INTEGER  N> 

PRINTER  I  1 6 
PRINT  PftGE 
PRINT  LINO) 

PRINT  •“  Name  the  record  which  w*11  receive  this  data." 

PRINT  LIN(1> 

PRINT  "  LENGTH  1  a ; N; "points. " 

F**" 

L INPUT  "  XXXXXX : F  or  XXXXXX 1 T 1 4 ",  Ft 

if  <F*C8;  no“F">  rind  <f*cs;  i  :<  -  >"T*  :•  *hen  2240 
INTEGER  Z(l:N) 

MftT  2*2  in 

IF  F ♦ £ §7 1 3 <  > ” F "  THEN  2320 

FCRE8TE  F I ,  I i)T  •:  2»N ' £58  1  ♦  INT<2*N  •  65538 )  +  3 
FFRINT  FS, 2<*  > 

GOTO  t'360 

!  DftTft  f  i  I  e : 

CREfiTE  FI,H  4 
ASSIGN  *1  TO  Ft 
PRINT  Si; 2  *> 

SUBEND 
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>  •  I  I  I  I  I  I  I  I  I  t  I  I  I  I  I  I  I  I  I  I  I  I  I  "OGRAM  "8-1 rt-l "  tllltlllllllttlllllllllllllllimil 

10  INTEGER  ZO:  30880), Z  rau< 1 : 38880) , Otservatt on*, N1 nt 

28  DIM  Raw  in*<8)£12J,ReiuU8<8)£l21 

21  (  <<<<<<<<<<<<<<<<<<<<<  PROMPTS  FOR  OPERATOR  »>>>>><>>>>>>>>>>>>>>>>>>>>>>>> 

38  PRINTER  IS  16 
48  PRINT  PACE 

90  PRINT  "  PROGRAM  8-in- l  removes  aircraft  Motion  froa  a  laser  prof i la" 

£0  PRINT  *  by  f i I  ter ing  the  data  in  thro#  at ops.* 

78  PRINT  LIN< 1 > 

80  PRINT  *  YOU  MILL  BE  REQUESTED  TO  PROVIDE  UP  TO  8  PAIRS  OF  FILES,* 

98  PRINT  *  EACH  CONSISTING  OF  A  RAM  DATA  FILE  FOLLOWED  BY  A  RESULT  FILE.* 
108  PRINT  LINO) 

118  PRINT  "  PLEASE  HAKE  FILE  NAMES  EXACTLY  6  CHARACTERS  LONG.* 

128  PRINT  LINO) 

138  PRINT  LINO) 

148  PRINT  "  TO  KEEP  TRACK  OF  THE  FILES  REQUESTED,  “ 

158  PRINT  *  ARRANGE  AND  Nt-IBER  YOUR  LIST  OF  FILE  PAIRS  BEFORE  CONTINUING.* 
1*0  PRINT  LINO) 

170  PRINT  *  Proas  CONT  to  go  or..* 

180  PAUSE 

198  FOR  Input -1  TO  8 
288  PRINT  PAGE 

218  PRINT  LINO) 

220  PRINT  “NAME  YOUR  <NEX7>  PAIR  OF  FILES.  <. DEVICE  CODE  MUST  BE  "tCHR 

a<34>t*F"tCMRS<34)t"  ).“ 

238  PRINT  LINO) 

248  PRINT  “  Exa*r-<*:  RawDaOF  CONT  Rosult:F8,l 

C0NT“ 

258  PRINT  LINO) 

268  PRINT  *  <  WHEN  FINISHED  NAMING  FILES,  ENTER  NOTHING  AND  PRESS  CO 

NT.  )" 

278  PRINT  LINO) 

288  PRINT  “FILE  PAIR  Input 

298  LINPUT  "  Raw  data  fi>o",Rau  inSOnput) 

388  LINPUT  "  Result  f  i  I  o" ,  RosuTt  SOnput  ) 

318  IF  (Raw  tn*< Input )>"* )  AND  (Roaul  tSOnput  >«*“  y  THEN  358 

320  IF  <Rau""i  n*<  Input  ■  £8}  1  X  >“F"  >  OR  (Rosult  SOnput  >C8;llO"F">  THEN  2 

98 

338  NEXT  Input 

348  !  !!!!!!I!!!!!!I!!!!!!!!!M!!>!!!!I!!!!'!!!I!!!!!!II!!!!!!!!I!II!!II!II!I! 

350  Jmaxolnput-l 


368  PRINT  PAGE 

376  PRINT  LINO) 

388  PRINT  "  YOU  WILL  BE  PROCESSING" ; J*ax; "FILE<S> .  " 

398  PRINT  LINO ) 

488  PRINT  "  Your  file<s>  are:  RAH  DATA  RESULT" 

410  PRINT  LINO) 

420  FOR  Q-l  TO  Jnax 

430  PRINT  "  "LPaw  in*<Q)t"  "tCHRS<240)t“  "iResul 1 *<Q ) 

448  NEXT  Q 

458  PRINT  LINO) 

468  PRINT  "  INSERT  NEW  DISC  IF  NECESSARY." 

478  PRINT  "  then  presi  CONT" 

488  PAUSE 

498  CALL  Warning  !  CRT  warning:  "DO  NOT  DISTURB  KEYBOARD” 

508  FOP  J* l  TO  Jn.ax  i  i  ••  i  .> ii ; t i !!!!!! M '••!!!!!!  i  !!!!!•!!•!!!!!!!!  H  ! 

510  Rau*“Rav  in*<J> 


50  <  <<<<<<<<<<<<<  READ  IN  DAT*  WITH  DISCONS  REMOVED  >>>>>>>>>)>>>>>  •>>.’>>>>.>>>> 


520 

FREAD 

Raws, Z_r aw<*  > 

538 

MAT  r 

■Z  raw 

540 

Obieruat  tor.a*ROU<Z) 

558 

Int  o*»R>:  j  ..  i  t  * 1  J ' 

551  »  «.:;«<«  CFEHTE  A  FILE  '  3  ACEPT  T  He  FINAL  PROFILE  >>>>>)>>>.>>>)  '■  /  > 

560  FCREATE  Into*,  INTOuS08*i  256  •♦InT-  30*00*2  65536 '*3 

578 

5.’1  '  <<<<<<<<V<  FEED  THE  RAW  I ATA  T-?'J  ThE  FIRST  ACTIVE  FILTER  >>>>  >/>»>>>)>> 
572  '  <<<<<<<<<<  <  REMOVE  'HE  HICH-FREOUEMCY  NOISE  J  >>>'>>>>>>>>>> 
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CALL  Act i ve_f i 1 ter_l <Z< O.Observat i on*>  !  LOUPASS  3  Hz 


t 


f 


HAT  Z  raw*Z  I  SINCE  HIGH-FREQUENCY  NOISE  HRS  BEEN  REMOVED 
<<<<<<<<<<  FEED  THE  RRU  DRTR  THRU  THE  SECOND  RCTIVE  FILTER  >>>>>>>>>>>>>>> 
<<<<<<<<<<  <  RERDY  DRTR  FOR  LOCRTINC  LOU  POINTS  )  >>>>>>>>>>>>>>> 


I 


CRLL  Ret  I ve_fi 1 ter_2<Z<«> , Obi ervat i on*>  !  HICHPRSS  .03  Hz 


I 


!  <<<< 


980 
590 
800 
SOI 
602 
610 
620 
630 
640 
650 
660 
662 
670 
680 
682 
690 
700 
710 
720 

730 

731 
740 
750 
760 
770 
780 

790  SUB  Low_poi  nt  <  INTEGER  Z_r»i.i<*> ,  INTEGER  2<* ) ,  Observat  i  on*,  lnt_l  engt h> 

800  INTEGER  2m i n< 1 : 1000) , Xmi n< 1 : 1000) , Interval , R1 ou , fthi gh, fl, 21 ou, No_of_i nt *, Bst 
art , Bit  op, B, 2_d> mi ni *h< 1 : 375) , C 


PRINT  PRCE 
PRINT  LIN<8) 

PRINT  "  LOW-POINT  ROUTINE  IN  PROGRESS." 

<<<<<  LOCATE  LOU  POINT!  AND  DRHU  LOW-POINT  ENVELOPE  >>>>>>>>>>>>>>>)>> 
CRLL  Low  point<2  raw< * > , 2< ♦  > , Ob*ervat i or*, 50) 

CALL  Warning  “  ”  !  CRT  warning:  "DO  NOT  DISTURB  KEYBOARD” 

<<<<<<<<<<<<<<<<  SMOOTH  THE  LOU-POINT  ENVELOPE  >>>>>>>>>>>>>>>>>>>>>>>>>>> 
CALL  Ret  i  ve_f  i  1  ter_3<2'*>  ,  Observat  i  on*)  *  LOUPASS  .13  Hz 


MRT  2*2  raw-2 
FPRINT  IntoF,  2<*) 


t 


Remove  aircraft  motion 
Store  final  profile 


I 


<<<■'<<<<<<<<<<<<  PROVIDE  A  HARD  COPY  OF  FINAL  PROFILE  >>>>>>>>>>>>>>>>>>>> 

CRLL  Graph i c*<Raui,  Into*,  Z<*> ,  Observat  i  ons) 


NEXT  J 


!!  |  j  |  !!)!!!!!!!  i  !  i)  !)!  )  MM)!  m  j  mm  I!)  j  !!!  |  !  M  )!!!!!!  I 


i 


END 


810 

820 

830 

8-*0 

8L0 

360 

870 

380 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

I860 

1070 

1030 

1090 

1100 

1110 

1120 

1130 


REAL  Del t a_x, Del ta  z 
MAT  Z*in*<0> 

Interval *0 
Flag-0 

Interval “Interval +1 

fll ow-< Interval -1 >*Int_l eng th+l 

Ahi gh* Interval *lnt_l ength 

IF  fihi gh<Obscrvat i on*  THEN  910 

F l ag* 1 

flhi gh*Observat i on* 

FOR  A*fllo‘w  TO  Rhigh 

IF  R-Alow  THEN  Zlow-C 
IF  2<R> >*21 ou  THEN  960 
21ou*2<fl) 

Xmi n< Interval )*A 
NEXT  A 

IF  ZlowOO  THEN  101O 

Zai n< Interval )*2_raw< Ahi gh >-Z<Rhi gh) 

X»in< Int  erval ;“Ahi  gh 
GOTO  1090 

IF  <Xmin<  Interval  >-10<- 1  >  OR  <Xmi n< Int erval >+10>Observat i on*)  THEN  1090 
Zl©w*Z_raw<Xmirr'Interval  >> 

FOR  C*Xai n< Interval )-10  TO  Xmi n< Int erval >♦ 10 


IF  2 


raw*  C>  >*21  op  THEN  1070 
Zloi.i»2_rau<C> 
Xmin< Interval >*C 


NEXT  C 

2*in<Interval . 
8  THEN  850 
int»*Interval 


=  3  raw  X. ■■  in'.  2  tv  erval  >  > 


IF  Flag 
No  of_ 

Flag-8 
Interval *8 
Interval “Interval ♦! 


1140 

IF  Int 

erv  a 

lot  THEN 

1170 

U50 

Bttar 

t*l 

1160 

GOTO 

1130 

1170 

IF 

Interval 

<  “No  _ 

1180 

IF 

Interval 

<■  No  C 

1 

1190 

IF 

Interval 

“No"* 

t  f 

1200 

IF  Bi 

top* 

3b*e  rvat i 

e— 

•H 

c 

o 

lEH 

1210 

i 

i  r,» 


»1  r-*E-l  B»t  art  “Xmi  n'  Interval -1 
♦  * ♦  1  THEN  Bit op*Xmt n< Int erval i 
ti»l  THEN  B*top*Obi«rvat ion* 

F ! ag* 1 
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1220 
1230 
1240 
1230 
1260 
1270 
1280 
1290 
1 3W0 
1310 
1320 
1330 
1340 
1330 
1360 
1370 
1330 
1390 
1400 
1410 
1420 
1430 
1440 
1430 
1460 
1470 
1480 
1490 
1300 
1310 
1320 
1530 
1340 
1350 
IS  50 
1570 
1330 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1730 
1760 
1770 
1780 
1790 
1800 
1810 
1620 
1830 
1840 
1850 
I860 
1870 
I860 
1890 
1900 
1910 


IF  Int*rval-1  THEN  0*1 1  *_z«Z» i  n<  l  t-Zjr  *u<  1 ) 

{  ~  ” 

IF  <Int*rv*1>l>  AND  <F1»g-0>  THEN  Dt )  t«_z-Z»>  n<  Int*rv*l  )-Z»in(  Interval -1  .> 
!  ™ 

IF  Ft*g»l  THEN  D*1 1  a_z-*._r au>:0b**rv*t  i ons>-Zni n(  Interval -1  > 

i  "" 

D*1 t*_x-B»top-B*t*rt  *1 
FOR  B«B*t*rt  TO  Bstop 

IF  Int#rv*l«l  THEN  Z<BT*Z_r *w< 1 >  +  D* 1 1 *_zxD*l t *_x*<B-B*t «rt -tl > 

IF  Interval >1  THEN  Z<B>*2iti  n<  Interval -1 >4-0*1 ta  z^Dvlta  x»<B-B*tart *1 > 

NEXT  B 

IF  F 1 ag-0  THEN  1130 
SUBEND 

SUB  Graphie*<Q*,P*, INTEGER  Z<*),0b*> 

INTEGER  A, B 
PLOTTER  IS  "GRAPHICS" 

GRAPHICS 

LINE  TYPE  3 

FOR  B-10  TO  60  STEP  10 
MOVE  0,B 
PLOT  0, B 
PLOT  120, B 
NEXT  B 

LINE  TYPE  1 
MOVE  0,85 

LABEL  "  Final  profit*  “4:P*fc“  <  f  row  "ibtt")  “ ;  Ob*;  “poi  nt*.  “ 

FOR  A=1  TO  Obs  STEF  10 
U*. 004*A 
V«Z<A>x33. 65+10 
IF  A*1  THEN  HOVE  U,V 
PLOT  U,V 

NEXT  A 

DUMP  GRAPHICS 
PRINTER  IS  0 

PRINT  LIN(S)  !  FEED  SOME  PAPER  WITH  PLOT 

PRINTER  IS  16 

SUBEND 

SUB  Act »ve_fi 1 t«r_l < INTEGER  2<  *> ,  ?v*erv*t  i  on*  > 

INTEGER  Z  *tor«<l:66>,St*tui,A,B,C,D 
OUTPUT  7,l;~AC!“ 

!  PREPARE  STEADY  VOLTAGE  FOR  INPUT 
FOR  A-l  TO  1000 

OUTPUT  7,  l;  “AB,2,  1,  1,  "«.VAl*<Z«'1  >•>*.“  !  " 

NEXT  A 

! 

!  INPUT  PROFILE 

i 

FOR  B-l  TO  Observation* 

OUTPUT  7,  l;  “AB,2.  1,  1 ,  "S  VAI.*<2 .  B)  >«-* ;  AI ,  t ,  1 , 1  !  " 

ENTER  7,1  BFHS  2  NOFORMAT ; St  u* 

ENTER  7,1  BFHS  2  NOFORMAT; Z> I  • 

NEXT  B 

!  CAPTURE  POINTS  DUE  TO  PHASE  LAG 
FOR  C-l  TO  66 

OUTPUT  7,  l;  "AB,  2,  1 ,  1 .  " t VAL*'  Z1. Obi* ru»t  i  on* >  >t“  ;  AI  ,1,1,1!“ 

ENTER  7,1  BFHS  2  NOFCRMAT; St  ttus 
ENTER  7,1  BFHS  2  N0FCPMAT;2  store- O 
NEXT  C 
SUBEND 

SUB  Act  iv*_fi  l  t*r_2< .INTEGER  2'  Obf  «>■-»*  i  on*  t 
INTEGER  2  »t  or  « *  1  \  6  .  S*  a*  -i  • .  A,  E,  C .  D 
OUTPUT  7,  l;~AC  " 

!  PREPARE  STEADY  VOLTAGE  FoP  INPUT 
FOR  A-l  TO  10OO 

OUTPUT  7, 1 ; "AB.2,2, 1, VALF-  Z'  I '  - '  " 

NEXT  A 

t  ' 

!  INPUT  PROFILE 

! 

FOR  B-l  TO  Observation* 
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i926 

193C 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2860 

2070 

2080 

2090 


OUTPUT  7, 1  {  “AB,  2 1 2,  1 ,  “iVALKZ'  BJ  )fc“;AI,l,2,  1  !  “ 
ENTER  7,1  BFHS  2  NOFORMAT ; St  at u* 

ENTER  7,1  BFHS  2  NOFORNAT; 2> B * 

NEXT  B 

i  CAPTURE  POINTS  DUE  TO  PHASE  LAG 

FOR  C-l  TO  66  _  ,  „ 

OUTPUT  7, 1 ;  "AB,  2, 2,  l ,  VAL*tZ<  Ob*  fru  at  i  on*)  )&“;«!, 1,2, 
ENTER  7,1  BFHS  2  NOFCF MAT;  St  %t  ..s 
ENTER  7,1  BFHS  2  NOFOFIAT; Z_st ore- C> 

NEXT  C 
SUBEND 

SUB  Act i ve_f i 1 t*r_3< INTEGER.  Z-  *> , Obser v at ions> 

INTEGER  Z_*Tor  e  <.  1 : 66  > ,  St  at  us .  A,  B ,  C ,  D 
OUTPUT  7, 1 ; " AC ! “ 

I  PREPARE  STEADY  VOLTAGE  FOR  INPUT 
FOR  A=1  TO  1000 

OUTPUT  7, 1 1  "AB,  2,3,  1 ,  “«YALF<Zf  1  >>&*'  1  " 

NEXT  A 


1  ! 


2100  > 

2110  !  INPUT  PROFILE 

2120  > 

?130  FOR  B»1  TO  Obser u*t  >om 

:•  1 40  OUTPUT  7,  l;  "AB,£.  3,  1,  “«.VAL*t  B)  )fc"  5  AI,  1 , 3, 1  '■  " 

"j50  ENTER  7,1  BFHS  2  NOFORMAT  1  St  at  u* 

2160  ENTER  7,1  BFHS  2  NOFOkMAi ; Z<B> 

2170  NEXT  B 

2180  !  CAPTURE  POINTS  DUE  TO  PHASE  LAG 

2190  FOR  C-l  TO  66  .  , 

?200  OUTPUT  7,  l;  "AB,2,3,  l,"iVAL*(2- :  Ot>  =  e  r '•  it :■  ions)  >«.“!  AI,  1,3,  1 

r?!0  ENTER  7,1  BFHS  2  NOFORMAT ;  St  at  us 

2220  ENTER  7,1  BFHS  2  NOFOF  'AT;  Z_*  >  ore  C  • 

2230  NEXT  C 

2240  FOR  D=1  TO  Observations 

'2*i0  IF  DC-Observat i ons-66  TmEN  2t,D>*2*. D*66!> 

2260  IF  D>Observat i ons-66  THEN  Z<D>=2_st ore < D- CObs er vat i ons-t 

2276  NEXT  D 

2280  SUBEND 

2290  SUB  Warning 

2500  PRINT  PAGE 

”310  PRINT  "  Aircraft  motion  being  removed. 

2320  PRINT  L!M(1) 

2330  PRINT  “ 

2340  PRINT  “ 

2350  PRINT  "  l 

2360  PRINT  " 

2370  PRINT  “ 

2380  FRINT 
2390  PRINT  “ 

2400  PRINT  " 

2410  PRINT  “ 

2420  PRINT  “ 

2430  PRINT  M 
i 440  PRINT 
2450  PRINT  “ 


2460  PPINT  " 

tl 

2470  PRINT  " 


2 460  PRINT  " 
2  4  ?0  PRINT  " 


2 ‘00  SOI END 
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